Quenching of phonon-induced processes in quantum dots due to electron-hole 

asymmetries 
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Differences in the confinement of electrons and holes in quantum dots are shown to profoundly 
impact the magnitude of scattering with acoustic phonons in materials where crystal deformation 
shifts the conduction and valence band in the same direction. Using an extensive model that 
includes the non-Markovian nature of the phonon reservoir, we show how the effect may be addressed 
by photoluminescence excitation spectroscopy of a single quantum dot. We also investigate the 
implications for cavity QED, i.e. a coupled quantum dot-cavity system, and demonstrate that the 
phonon scattering may be strongly quenched. The quenching is explained by a balancing between the 
deformation potential interaction strengths and the carrier confinement and depends on the quantum 
dot shape. Numerical examples suggest a route towards engineering the phonon scattering. 



It is well-known that cavity QED systems undergo de- 
phasing due to interactions with the environment, and 
that this leads to loss of quantum mechanical coherence 
between the different states of the system. In particu- 
lar, solid-state based cavity QED systems such as semi- 
conductor micropillars [lH3| and photonic crystal cavi- 
ties @, @ are strongly affected by dephasing induced by 



phonon scattering [61-113] . Phonon-induced decoherence 



thus makes the observation of such effects as vacuum 
Rabi oscillations [lij much more difficult than in atomic 
cavity QED systems and impairs the realization of a scal- 
able solid-state platform for quantum information tech- 
nology. For instance, the requirement of indistinguisha- 
bility of subsequent emission events from a single-photon 
source places strict limitations on the amount of dephas- 
ing that can be accepted [15l4l7|. 

In this Letter we investigate the role of carrier con- 
finement on phonon scattering in a semiconductor quan- 
tum dot (QD). We employ a comprehensive theoretical 
model that takes into account the non-Markovian na- 
ture of the phonon reservoir [l2j , thus avoiding the stan- 
dard approach of describing phonon-induced decoherence 



by a pure dephasing rate 18 20]. Surprisingly, we find 



that acoustic deformation potential scattering may be 
completely quenched under certain conditions, depend- 
ing on the degree of confinement of the involved elec- 
tronic states. The effect is investigated in detail for two 
cases. 

Firstly, we consider the case of photoluminescence 
spectroscopy of a single QD and find that the lumines- 
cence is suppressed at certain detunings due to a quench- 
ing of phonon scattering. For spherical confinement po- 
tentials we derive approximate analytical results for the 
scattering rate. These results show that the effect orig- 
inates from the difference in the spatial confinement of 
electrons and holes, usually neglected in theoretical treat- 
ments. 

Secondly, we consider a coupled QD-cavity system. We 
find that the quenching of the phonon scattering strongly 



affects pure dephasing with consequences e.g. for the 
indistinguishability of single photon sources. We further 
expand the description to realistic QD structures, which 
are analyzed numerically, and the conditions for reducing 
phonon scattering are established. 

Other approaches towards controlling the degree of 
phonon scattering include the use of phononic bandgap 
structures for suppressing vibrational modes [2l| . Struc- 
tures with simultaneous photonic and phononic bandgaps 
have been discussed theoretically [22j and experimentally 
realized [23[ with phononic bandgaps in the GHz-regime. 
In cavity QED, the most relevant phonons belong to the 
acoustic branch and have energies in the THz-regime, 
hence the present phononic band gap structures are pri- 
marily of interest for improving the optomechanical cou- 
pling [24]. Another approach towards manipulating the 
electron-phonon interaction has been demonstrated by 
placing the QD near a surface, thereby changing the 
phononic dispersion [jjj ]. 

At first we consider a two-level QD with transition en- 
ergy hujQD illuminated by a CW laser with frequency 
wl in a standard photoluminescence excitation experi- 
ment. We limit ourselves to detunings below 10 meV, 
where the the deformation potential coupling between 
electrons and longitudinal acoustic phonons constitutes 
the dominating coupling mechanism [26j. Initially we 
consider spherical-parabolic confinement potentials, giv- 
ing ground state wavefunctions 



<M r ) = T 



exp [-r 2 /(2zO] , (1) 



with v £ {e, h} representing the electron and the hole re- 
spectively. Numerical results for more realistic QD struc- 
tures will be presented later. The differences in carrier 
confinement lead to wavefunctions with different effec- 
tive widths for electrons (w e ) and holes (wj,). The equa- 
tions of motion for the system are derived using a quan- 
tum master equation approach [27| exploiting recent re- 
sults [IH for the reduced density matrix, p(t), where the 
phononic degrees of freedom have been traced out. The 
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FIG. 1: (color online), (a) Stationary population of a bare 
QD, when excited by a CW field, corresponding to a Rabi- 
frequency Ml = 10 fieV. The black curve describes the case 
without phonon-assisted coupling, and the QD decay rate is 
r = Ins -1 , (b) The effective phonon spectrum for spherical 
wavefunctions for varying wh and w e . For both plots, + w\ 
is kept constant, and T = 4K. 

system evolves according to 

d tP {t) = -m- l [H s ,p{t)]+S hA [t) + L{t), (2) 

where H s — h(ujQ^~Lu-^)a ee + hfl(a eg +a ge ) is the system 
Hamiltonian for the QD and the optical excitation [35j |. 
corresponding to a Rabi- frequency Ml = 10 /ieV. Sla(*) 
is a time-local phonon- induced scattering term, and L(t) 
represents losses included through the Lindblad formal- 
ism 



28] 



The latter accounts, through the rate V, for the 
decay of the excited QD state due to spontaneous radi- 
ation into other modes than the cavity mode as well as 
non-radiative processes. 

Solving Eq. @ numerically Hi}, we calculate the sta- 
tionary population of an InAs QD embedded in GaAs 
in depending on the detuning between the laser and the 
QD ground state resonance, see Fig. [T]i. Comparing 
the black and blue curve in Fi g. QJi shows that the 
QD indeed, as is well-known (29l. 130} . may be populated 
by a phonon-assisted process, even when the laser and 
the QD are off-resonant. However, when the widths of 
the electron and hole wavefunction differ, as for the red 
curve, the QD population and thus the photolumines- 
cence intensity, is suppressed at a specific detuning. As 
we shall show, this implies that the phonon scattering 
is quenched, which introduces a new and interesting ap- 
proach for manipulating the electron-phonon interaction. 

The strength of the carrier-phonon scattering is quan- 
tified by the effective phonon spectrum, T)(u), describ- 
ing how the phonon modes interact with the QD at a 
given temperature. Positive phonon frequencies refer to 
a phonon of energy fiw being emitted into the surround- 
ings, whereas negative phonon frequencies indicates the 
absorption of a phonon of energy h\uj\ from the surround- 
ing environment. T>{ui) is calculated as the real part of 
the Fourier transform of the phonon reservoir correlation 
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n k <5(u; + Wk) + ("k + 1)8(uj - Wk) 



(3) 

where k is the phonon wavevector, = 
l/[exp(fiwk/A;BT) — 1] is the average thermal occu- 
pation number of the phonon mode k at temperature T. 
Af k = M^ e — M^ h is the effective electron/hole-phonon 
coupling matrix element with 



The electronic form factor, 



hk 

2dc t V 



•7v(k). 



TvQl) = / dr|0„(r 



l 2 e ikr , 



(4) 



(5) 



describes how a phonon with wavevector k interacts with 
the confined carriers, governed by the overlap between 
the carrier and phonon wavefunctions. We assume bulk 
phonons with a linear dispersion relation, lu^ — q|k|. D e 
(Dh) is the deformation potential constant of a conduc- 
tion band electron (valence band hole), c; is the velocity 
of longitudinal acoustic waves, d is the bulk density, and 
V is the phonon quantization volume. 

For spherical wavefunctions, cf. Eq. (JTJ, the effective 
phonon spectrum reduces to 



Airdc)" 1 - e-P huj 



(6) 



and is shown in Fig. QJ>. The magnitude of T>(uS) is 
generally smaller for to < compared to u> > due to 
the low probability of thermally excited phonons. From 
Eq. (|5|) we note that a spatially narrow wavefuntion is 
wide in k-space and thus interacts with many phonon 
modes, explaining why T>(uS) broadens and increases in 
magnitude as wh decreases. For T > 0, T>(lu) has three 
zeros at 



0, u 2 



Acf 



"o In ! 
• 2 * D h 



(7) 



In materials where D e / > 0, like GaAs [g, [l8|, dips 
thus appear in the effective phonon spectrum for non- 
zero phonon energies for w e ^ w^. The zeros in the 
effective phonon spectrum appear exactly when |M k | oc 
\D e T e {k = u/ci) - D h T h {k = uj/ci)\ = => D e T e {k = 
u/c{) = D h T h {k = uj/ci). 

Based on this we can explain the quenching of phonon- 
induced processes. The deformation potential interaction 
occurs in general due to the different values of the defor- 
mation potential constant in the conduction and valence 
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FIG. 2: (color online). Coupled cavity-QD system interact- 
ing with a phonon reservoir. The optical coupling strength is 
g, the QD-cavity detuning is A (positive for the case shown), 
and M ee /h>i are electron/hole-phonon coupling matrix ele- 
ments. F is the QD population decay rate, and k is the leak- 
age rate from the optical cavity. 



band. However, for a specific detuning, this can be com- 
pensated by a difference in the confinement of the elec- 
tron and the hole through the form factor in Eq. (|5|). 
fn the case of identical electron and hole wavefunctions, 
w e = Whi the zeros appear at infinity, and the effect is 
not apparent. 

The quenching of phonon processes has interesting con- 
sequences for cavity QED. Instead of illuminating the QD 
with a laser, the QD is now placed inside a single-mode 
optical cavity, detuned by A = wqd — w cav from the cav- 
ity resonance, see Fig. [5] In this case, H s in Eq. @ 
is replaced by the combined QD-cavity system Hamilto- 
nian (corresponding to replacing 57 with the optical in- 
teraction strength g with the cavity mode), see e.g. [lij . 
Furthermore, the escape of photons through the cavity 
is included via a rate k in the Lindblad term L(t) in Eq. 

By solving Eq. ([2]), we obtain the QD decay curves 
shown in Fig. [3^- The structure of the phonon bath 
results in an asymmetry, expressed in the possibility of 
the QD to couple to a red-tuned (lower energy, A > 0) 
cavity by the emission of an acoustic phonon, but a lack 
of coupling to a blue-tuned (higher energy, A < 0) cavity 
by phonon absorption at low temperatures, where the 
population of thermally excited phonons is low [(| fl2| . 
This is in accordance with T>(ui) in Fig. [T}d. In the limit 
of |A| — > oo, the QD and the cavity decouple, and the 
QD decays with the rate T. By single-exponential fits 
to the QD decay curves, the lifetime of the excited QD 
state is calculated, see Fig. [3Jd. From this we determine 
the lifetime ratio ta<o/ta>o, see Fig. [3fc. Surprisingly, 
we observe that when w e ^ Wh, a non-zero QD-cavity 
detuning exists, where no lifetime asymmetry is present. 

To examine the physical origin of the lifetime asymme- 
try, we consider the total decay rate of the excited QD 
state in the large detuning limit A 3> g 15j, 



r tot «r + 2. 9 2 ^i- 



7t 2 ot + A 2 



t 



-V{oj = A) 



(8) 
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FIG. 3: (color online), (a) Time evolution of the population 
of an initially excited QD for different detunings. The solid 
(dashed) curve is for negative (positive) detuning, and the 
electron and the hole are equally confined, w e = Wh = 5nm. 
The temperature is T — 4K, and we use fig = 150 fieV, T = 
f ns _1 , and hn = fOO /ieV. (b) Lifetimes, r, of the excited QD 
state plotted versus detuning for different values of w e and u>h , 
keeping the volume parameter w'i + w\ constant. The solid 
(dashed) curves indicate negative (positive) detunings. (c) 
The degree of asymmetry versus detuning. 



has two contributions besides the decay rate T. The first 
term in the square brackets represents the usual Purcell 
enhancement rate, while the second term represents the 
decay of the QD through the cavity by the simultaneous 
emission/ absorption of a phonon. Thus the behaviour of 
T>{uS) translates directly into the behaviour of the life- 
time ratio plot in Fig. [3J:, such that non-zero QD-cavity 
detunings exist, where phonons do not affect the lifetime 
of the excited QD state. 

The suppression of the effective phonon density is ex- 
pected to affect not only the QD lifetime, but all phonon- 
induced effects such as pure dephasing. The degree of 
single-photon indistinguishability [f6| is typically quan- 
tified as the ratio between the coherence time and life- 
time of the emitter [l7j]. For our system, this ratio is 
r t ot/(r t ot + 27 ph ), with 7 ph being the pure dephasing 
contribution from the phonons. This hints at the im- 
portance of having a high decay rate compared to the 
dephasing rate. 

The phonon-induced pure dephasing rate [33] is calcu- 
lated from the phonon bath correlation function, 

7 P h(i) = Rc j/r 2 (l -K)J dt' V(t') 

+h- 2 K J dt' cos (tVA a +4$ a ) P(*')}, (9) 
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Wh — 5 nm 
Wft = 4 nm 
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FIG. 4: (color online), (a) Time evolution of phonon- induced 
dephasing rate for different QD-cavity detunings. Inset show 
the long-time value of the dephasing rate versus detunings, 
Wh = 3nm, w e = 6nm, and T = 4K. (b) Time evolution 
of phonon-induced dephasing rate for different wavefunction 
widths, keeping w^ + w^ constant, and A = lmeV. The inset 
quantifies the short-time scattering by showing the integrated 
value of 7 P h(t) from t — to 3 ps for different wavefunction 
widths. 



where if is a parameter depending solely on the ratio 
A/g. For bulk phonons, T>(t) decays within 5 ps fl5j j . 
As seen from Fig. 2^, the initial behaviour of 7 P h(i) is 
governed by the bare QD-phonon coupling, whereas the 
long-time value depends on the cavity detuning. On the 
short time scale, the phonons participate in non-energy 
conserving virtual processes. Thus the electrons interact 
with the full phonon spectrum D(lu). The long-time (t — > 
oo ) value of 7 P h(i), on the other hand, implies a Fourier 
transform of V(t) and depends only on T>(u = 0) and 
T>(cj = ±-\/A 2 + 4g 2 ), where the former is zero. 

The phonon-induced dephasing of the QD-cavity sys- 
tem may be reduced by minimizing both the short-time 
and long-time scattering. The former may be achieved 
by reducing the overall amplitude of T>(uj), as illustrated 
in Fig. [4Jd for Wh ~ 3.4 nm, compare with Fig. [TJd. The 
reduction of long-time scattering is obtained by using a 
QD-cavity detuning such that y/ A 2 + Ag 2 coincides with 
a dip in T>(lu), as seen from the inset in Fig. It turns 
out that for the indistinguishability of single photons, the 
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FIG. 5: (color online), (a) The effective phonon spectrum cal- 
culated from the FEM wavefunctions of the three QDs shown 
in the insets with matching colors, T = 4K. All QDs have 
constant volume and side wall slope, (b) The effective phonon 
spectrum for the shallow dot in (a) for different material com- 
positions, In^Gai-yAs, of the QD. 



reduction of the overall amplitude is the most important 
& 

QD structures realized by epitaxial growth are not 
spherically symmetric [H} , and to investigate the quench- 
ing for more realistic structures, we model truncated con- 
ical QD structures by solving the one-band effective mass 
Schrodinger equation using the Finite Element Method 

The effective phonon spectrum for the FEM wavefunc- 
tions is shown in Fig. [5^ for an InAs QD and wetting 
layer embedded in a GaAs barrier 38[. The overall am- 



plitude of the spectrum is the largest, when the carriers 
are strongly confined in two directions, i.e. for the tallest 
dot, compared to the shallow dot where the carriers are 
only strongly confined in one direction. These results are 
qualitatively explained by the form factor in Eq. ([5]). 

To examine the role of material composition, we con- 
sider an InyGai_yAs QD, where the amount of gallium 
in the QD and wetting layer is varied. By increasing the 
gallium concentration the band offsets shrink 34| and the 
effective mass changes, which results in the phonon spec- 
trum in Fig. [SJd. Thus, w e and Wh both increase, result- 
ing in a lower overall amplitude of the effective phonon 
spectrum. Due to the heavier hole mass, the asymme- 
try between the electron and hole wavefunctions also in- 
creases, moving the dip in the spectrum towards lower 
frequencies, as predicted by the spherical model in Eq. 
([7]) when w 2 — w\ increases. The position of the dips 
may in this way be changed depending on the growth 
parameters of the QD. For small gallium concentrations, 
the confinement energies may become comparable to the 
Coulomb energy, in which case the exitonic nature of the 
electron-hole pair will have to be taken into account. 

In conclusion, through the analysis of a comprehensive 



5 



model, we have learned that electron-phonon scattering 
may be suppressed due to differences in the spatial con- 
finement of electrons and holes. We suggest that the 
effect may be measured by photoluminescence excitation 
spectroscopy, but also that it should strongly affect the 
decoherence properties of cavity QED systems. We ex- 
pect this phonon-quenching to be important for the en- 
gineering of the phononic properties of devices for quan- 
tum information technology, like single-photon sources 
and switches. 
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